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INTRODUCTION 

This  research  addressed  problems  of  combustion  instability  in  liquid-propellant 
rocket  motors,  with  attention  given  to  both  hydrogen-oxygen  systems  and  storables. 

The  studies  were  theoretical  and  involved  both  asymptotic  and  statistical  methods 
of  applied  mathematics.  The  aim  w'as  to  improve  understanding  of  the  phenomena 
of  combustion  instability  and  to  indicate  potentially  useful  approaches  to  instability 
analyses. 

OBJECTIVE.S 

The  overall  objective  of  this  research  was  to  improve  understanding  of  the  mecha¬ 
nisms  by  which  flow,  mixing  and  combustion  process  are  coupled  to  acoustic  fields  in 
liquid-propellant  rocket  motors.  Particular  attention  wcis  focused  on  effects  of  spatial 
and  temporal  inhomogeneities  of  the  acoustic  media  which  are  associated  with  tur¬ 
bulence  and  with  two-phase  flow.  Appropriate  statistical  approaches  were  provided 
for  different  types  of  inhomogeneities.  In  addition,  amplification  mechanisms  coupled 
with  finite-rate  chemical  reactions  were  analyzed  by  use  of  activation-energy  asymp¬ 
totics  and  other  eisymptotic  methods. 

ACCOMPLISHMENTS 

Rotational  Inviscid  Flow  Effects  in  Laterally  Burning  Motors 

As  a  first  step  towards  investigating  instabilities,  steady-state  flow  fields  in  rocket 
motors  were  addressed.  A  theoretical  analysis  to  determine  the  effects  of  mass  addi¬ 
tion  on  the  inviscid  but  rotational  and  compressible  flowfield  in  a  porous  duct  with  the 
injection  rate  dependent  of  the  local  pressure  was  performed  for  large  ratios  of  length 
to  duct  diameter.  The  problem  of  describing  the  flow  was  reduced  to  the  solution  of  a 
single  integral  equation.  The  ratio  of  specific  heats, 7,  and  a  constant  pressure  expo¬ 
nent,  n,  measuring  the  dependence  of  the  rate  of  mass  injection  on  the  local  pressure, 
are  the  parameters  of  the  solutions.  The  integral  equation  w’as  solved  numerically, 
and  parametric  results  were  obtained  for  7  varying  from  1  to  |  and  for  n  varying 
from  0  to  1.  A  choking  phenomenon  is  exhibited  at  a  critical  length  of  the  duct  in  the 
vicinity  of  which  the  Mach  number  approaches  unity.  The  choking  condition,  which 
is  relevant  to  the  operation  of  nozzieless  solid-propellant  rocket  motors,  was  obtained 
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parametrically  in  the  present  study  and  compared  with  corresponding  results  for  irro- 


tational,  quasi-one-dimensional  flow.  The  rotationality  reduces  the  choking  pressure. 


-‘Iity  Codes 


Effects  of  Spatial  and  Temporal  Inhomogeneities  Produced 
by  Turbulence  and  by  Two-Phase  Flow 
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Iwo  distinct  approacnes  to  tne  tlieory  oi  stocnastic  enects  in  tne  combustion  in- - 

stability  of  liquid-propellant  rockets  were  developed.  In  one  approach,  dealing  with 

effects  of  spatial  inhomogeneities  in  the  acoustic  medium  associated  with  two-phase 
flow  and  wdth  turbulence,  attention  is  restricted  to  the  particular  class  of  interactions 
in  which  the  characteristic  time  scde  of  turbulence  is  asymptotically  large  compared 
with  that  of  acoustic  waves,  so  that  the  acoustic  medium  is  considered  to  be  a  qua¬ 
sistationary  random  medium.  In  addition,  it  is  assumed  that  the  acoustic  wavelength 
is  long  compared  with  any  scale  of  inhomogeneity  within  the  chamber,  motivating 
a  homogenized  description,  in  which  averages  of  state  variable  becomes  appropriate. 
On  the  basis  of  these  assumptions,  the  Navier-Stokes  equations  for  two-phase  flows 
were  reduced  to  a  non  homogeneous  stocheistic  Helmholtz  equation.  Source  terms  of 
the  equation  were  identified  as  arising  mainly  from  phase  change,  from  homogeneous 
chemical  reaction,  and  from  spatial  variations  of  wave  properties  such  as  sound  speed. 

To  handle  the  stocheistic  wave  equation,  the  method  of  smooth  perturbation  was 
used  and  resulted  in  a  Helmholtz  equation  for  an  equivalent  deterministic  acoustic 
medium.  With  use  made  of  Green’s  integral  theorem,  the  dispersion  relation  for 
transverse  acoustic  modes  in  cylindrical  chambers  was  obtained  in  the  form  = 
^0  +  +  ^ob  •  k  -h  k^c,  where  a,  b  and  c  are  calculated  from  appropriate  averages  of 

source  terms.  In  particular,  the  k^c  term  corresponds  to  the  stochastic  contributions 
arising  mainly  from  the  spatial  correlations  of  the  sound  speed  with  the  pressure 
response  function,  and  of  the  local  Mach  number  with  the  velocity  response  function. 
Compared  with  the  deterministic  contributions  to  the  linear  growth  (or  damping) 
rate,  the  stochastic  contributions  are  found  to  be  of  the  order  of  where  £t 

is  the  characteristic  length  scale  for  turbulence  and  R  is  the  chamber  diameter,  thus 
implying  that  most  of  the  stochastic  contributions  comes  from  the  larger  turbulent 
eddies.  The  stochastic  effects  are  expected  to  be  significant  in  the  transition  regime 
in  which  the  deterministic  growth  rate  is  close  enough  to  zero  for  transition  between 
stable  and  unstable  domains  to  occur  reasonably  frequently  by  stochastic  variations 
of  the  growth  rate. 
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In  another  approach,  emphasizing  influences  of  turbulent-induced  noise,  the  be¬ 
havior  of  high-frequency  combustion  instabilities  influenced  by  random  temporal  vari¬ 
ation  of  the  linear  growth  rate  in  this  same  transition  regime  are  addressed.  The 
Rayleigh  criterion  was  generalized  to  account  for  turbulence-related  spatial  variations 
of  the  combustion  response.  Special  attention  was  given  to  the  distinguished  limit  in 
which  the  acoustic  oscillations  are  rapid  compared  with  the  turbulent  fluctuations, 
w’hich  in  turn  are  rapid  compared  with  the  linear  growth  rate  of  the  instability.  It  wa^ 
shown  that  nonlinear  acoustics  of  the  instabilities  need  to  be  considered  to  describe 
the  influences  of  turbulence.  When  this  nonlinearity  involves  supercritical  bifurcation, 
the  turbulence  tends  to  decrease  the  most  probable  intensity  of  the  instability.  When 
it  involves  subcritical  bifurcation,  the  turbulence  can  produce  a  bimodal  probability 
density  function  for  the  intensity  of  the  instability,  with  appreciable  probabilities  of 
high-amplitude  and  low-amplitude  acoustic  oscillations  but  small  probabilities  of  os¬ 
cillations  of  intermediate  amplitudes.  These  phenomena  can  have  a  bearing  on  erratic 
pressure-amplitude  bursts  sometimes  observed  in  liquid-propellant  engines. 

A  key  observation  of  this  last  study  is  that  the  turbulence  influences  on  the  acous¬ 
tics  appear  as  multiplicative  rather  than  additive  noise.  Comparison  of  the  turbulence 
coherence  time  with  the  characteristic  time  of  the  combustion  instabilities  identifies 
two  distinct  interaction  classes,  depending  on  whether  the  coherence  time  is  or  is 
not  much  smaller  than  the  characteristic  instability  growth  time.  Larger  instabili¬ 
ties  are  associated  with  the  second  of  these  classes  while  the  transition  behavior  in 
the  first  class  involves  random  walks  of  the  stochastic  linear  growth  rate,  requiring 
a  statistical  description.  In  a  white-noise  approximation,  the  Stratonovich  rather  Ito 
interpretation  of  the  stochastic  integral  must  be  employed,  and  the  Markov  property 
of  the  white-noise  process  then  enables  the  evolution  of  the  probability  density  of  the 
pressure  amplitude  to  be  described  by  a  Fokker- Planck  equation,  whose  stationary 
solutions  lead  to  transition  behavior  that  depends  on  the  character  of  the  bifurca¬ 
tion.  There  is  a  bimodal  behavior  for  subcritical  bifurcation  in  the  hysteresis  region. 
The  size  of  the  bimodal  transition  region  is  proportional  to  the  timo  integral  of  the 
autocorrelation  function  of  the  turbulence-induced  fluctuations  of  the  linear  growth 
rate.  This  result  provides  a  possible  mechanism  for  resurging  patterns  of  the  pressure 
amplitude. 
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Amplification  Mechanisms  Coupled  with  Finite-Rate  Chemistry 

Many  different  processes  can  contribute  to  linear  amplification  of  acoustic  oscil¬ 
lations  in  liquid-propellant  rockets,  and  asymptotic  analyses  of  these  processes  can 
help  in  assessing  combustion  instabilities.  It  is  the  intent  of  the  present  investiga¬ 
tion  to  analyze  a  number  of  these  processes.  One  such  process  is  the  strained  planar 
diffusion  flame  considered  here.  Previous  studies  of  diffusion-flame  response  postu¬ 
lated  infinite  chemical  reaction  rates  (the  Burke  Schumann  approximation),  while  the 
present  work  is  taking  into  account  the  influences  of  finite-rate  chemistry,  which  can 
become  important  in  environments  having  high  turbulence  intensities.  The  results 
demonstrate  that  the  high  sensitivity  of  the  chemical  reaction  rates  to  temperature 
fluctuations  can  underlie  important  amplification  mechanisms.  As  an  initial  simpli¬ 
fication,  a  one  step,  irreversible  Arrhenius-type  chemical  reaction  rate  is  emploj’ed, 
and  a  gciseous  counterflow  diffusion  flame  is  adopted  to  represent  flamelets  subjected 
to  nonuniform  flow  fields  caused  by  turbulent  fluctuations.  It  is  intended  later  to 
introduce  rate-ratio  asymptotics  for  hydrogen-oxygen  systems. 

The  analysis  was  performed  by  activation-energy  asymptotics.  The  resulting  flame 
structure  for  a  given  value  of  the  reaction-sheet  location  is  described  by  two  sets  of 
equations,  one  for  the  transport  of  momentum,  thermal  energy  and  reactant,  and  the 
other  for  the  corresponding  overall  reaction  rate  of  the  flame.  The  acoustic  response 
of  the  flame,  obtained  from  the  linear  analysis,  is  found  to  be  determined  by  two 
mechanisms,  namely,  oscillations  of  the  reaction  sheet  induced  by  acoustic-produced 
fluctuations  of  the  reaction  rate,  and  oscillations  of  the  field  variables  produced  by 
the  transport-zone  response.  The  results  show  that  analyses  for  the  acoustic  response 
of  flames  that  do  not  consider  finite  reaction  rates  could  significantly  underestimate 
the  amplification  rate.  Future  works  will  apply  similar  methods  to  different  types  of 
flamelets  which  may  be  more  realistic  for  rocket-engine  combustion,  such  as  single¬ 
droplet  combustion,  with  and  without  forced  convection  and  supercriticality. 

CONCLUSION 

This  research  has  contributed  to  our  understanding  of  ways  to  describe  combustion 
instabilities  in  liquid-propellant  rocket  motors. 
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